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SPECIFICATION 
TITLE 

"METHOD FOR SURFACE-CONTOURING OF A THREE-DIMENSIONAL 

IMAGE" 

BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention concerns a method for surface contouring of a 
three-dimensional image of a subject. 
Description of the Prior Art 

For various applications (for example photographic analysis, industrial 
image processing, monitoring, traffic systems, media technology and 
animation), image processing and subj e ct object recognition of two- 
dimensional, and to a certain extent three-dimensional, images of subj e cts 
objects have progressed relatively far. The image data sets forming the 
images generally are , as a ru le , relatively w e akly d i srupt e d undisturbed and 
are acquired, for example, with optical sensors (see, for example, Inspect. 
MessTecAutomation, special publication 10/01, GIT Publishing Darmstadt, 
2001, or Vision Systems Design, Vol. 6 #9, Pennwell, Sept. 2001). 

If the image data sets are relatively w e akly d i srupt e d undisturbed , the 
surface contours of the images of the subj e cts objects can be relatively simply 
determined by means of known methods. With linear or nonlinear filtering, 
noise in the image data sets can be further reduced. In order to recognize 
textures of the images, pattern matching can be implemented. Bord e rs of the 
i mages Object borders can be determined, for example, by 2D- or 3D gradient 
format i on calculation . Techniques based on gradient formation calculation 



are, for example, "active contours", specified in Chalana V., Linker D., "A 
Multiple Active Contour Model for Cardiac Boundary Detection on 
Echocardiographic Sequences", IEEE Trans. Med. Imaging 15, #3, June 
1996, or what ar e known so-called as "snakes", specified in Mclnerney T, 
Terzopoulos D: "T-snakes: Topology Adaptive Snakes", Medical Image 
Analysis 4, pages 73-91, 2000. Imaged subj e cts objects of the same 
brightness, color or texture can, on the other hand, be "filled out", for example 
by a technique known as region growing, and thus be completely described 
for a specific task (see Jendrysiak U., "Segmentierung von Schnittbildern". 
Spektrum der Wissenschaft, Dossier 1/1999: Perspektiven in der 
Medizintechnik 24-29, 1999). 

In medical technology, however, it may be the case that imaged 
healthy and diseased tissue exhibit a relatively small contrast difference, and 
image data sets are ro l at i vo l y sovoroly very noisy. This is true in particular for 
x-ray images (computed tomography, angiography, fluoroscopy), images 
produced with a magnetic resonance device or in nuclear medicine, and m 
part i cu l ar particularly for images based on ultrasound (compare Marais P., 
Brady J., "Detecting the Brain Surface in Sparse MRI Using Boundary 
Models", Medical Image Analysis 4, 283-302, 2000, Chalana V., Linker D., "A 
Multiple Active Contour Model for Cardiac Boundary Detection on 
Echocardiographic Sequences", IEEE Trans. Med. Imaging 15, #3, June 
1996, or Sakas G., "Dreidimensionale Bildrekonstruktion aus Ultraschall- 
Daten", Spektrum der Wissenschaft, Dossier 1/1999: Perspektiven in der 
Medizintechnik 18-24, 1999). 
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In part i cular i n a diagnost i c find i ng Particularly when searching for a 
diagnosis for a patient, fine subj e ct object differences must be differentiated in 
the overall depth of the body or of the organs of the patient. Since the image 
data sets acquired with medical devices often are s e v e r el y very noisy, and 
acquired structures can be relatively complex, a determination of a contour in 
a two-dimensional image is already relatively computation-intensive. The 
determination of a contour of a three-dimensional image of a subj e ct an object 
is even more difficult. Known methods therefore have been refined, and often 
ad hoc adapted ad hoc to[[,]] the present problem. An example is specified in 
Chalana V., Linker D., "A Multiple Active Contour Model for Cardiac Boundary 
Detection on Echocardiographic Sequences", IEEE Trans. Med. Imaging 15, 
#3, June 1996. Another relatively computation-intensive method is specified 
in Jendrysiak U., "Segmentierung von Schnittbildern". Spektrum der 
Wissenschaft, Dossier 1/1999: Perspektiven in der Medizintechnik 24-29, 
1999. 

SUMMARY OF THE INVENTION 

It is an object of the invention to provide a method with which a surface 
contouring of a three-dimensional image of a subject an object can be 
implemented with less calculation effort. 

This object is achieved in accordance with the invention by a method 
for surface contouring of a three-dimensional image of a subject including the 
steps of placing a first set of planes through the three-dimensional image, 
determining the contours that are imaged in each of the planes of the first set 
of planes and are associated with the surface contour of the three- 
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dimensional image, and m e rging composing the contours determined in each 
plane of the first set of planes into a surface grid associated with the surface 
contour of the three-dimensional image. 

In accordance with the invention the surface of a three-dimensional 
image of a subject shou l d is to be contoured or segmented. The three- 
dimensional image is stored in a volume data set, thus a data set is formed of 
a 3D b l ocks block B(x,y,x) with KxLxM voxels (volume elements). The three- 
dimensional image, alternatively can be stored in the form of a data set that 
contains K slices of images with LxM pixels (picture elements) per image. 

For this, the planes of the first set of planes are placed through the 
three-dimensional image. The contour (that is associated with the surface 
contour of the three-dimensional image) is subsequently determined for each 
one of these planes. The determination of the contour associated with the 
individual planes corresponds to the determination of contours of two- 
dimensional images. After this, the individual contours are merged composed 
into a surface grid associated with the surface contour of the three- 
dimensional image. 

In a preferred embodiment of the invention, the planes of the first set of 
planes intersect in a first straight line that, according to a particularly 
advantageous variant of the invention, proceeds substantially through the 
geometric center (center of grav i ty mass ) of the three-dimensional image. 
This is particularly advantageous givon for a convex three-dimensional image. 

If, for example, the volume data set forming the three-dimensional 
image exhibits a relatively weak contrast, the volume data set can thus 
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undergo a filtering, in particular for an automatic localization of the three- 
dimensional subj o ot object (according to a variant of the invention) for the 
determination of the geometric center of the three-dimensional image. The 
filtering in particular can be a median filtering, and the geometric center in 
particular can be automatically determined. The volume data set thus is 
significantly smoothed while preserving edges. Via such a significant 
smoothing, the detailing of the boundaries is may be partly lost, however the 
geometric center of the three-dimensional image can be r el at i v e ly precisely 
determined. 

If the subject to imaged is a part of a living organism, and the three- 
dimensional image is acquired by means of a medical technology device, the 
geometric center of the three-dimensional image eae may be determined 
during a navigation-guided treatment of the living organism, according to a 
further embodiment of the invention. The geometric center can be 
determined, for example, during the navigation-guided treatment by the 
position of the focus of a navigated medical instrument. 

A medical device suitable for the navigation is, among other things, 
specified in German OS 199 51 502. For example, a volume data set of the 
region of interest of the body of a patient is generated with an x-ray device. In 
particular during a minimally invasive procedure, for example, a surgeon 
inserts the medical instrument into the body of the patient. A position sensor 
of the navigation system is arranged at the medical instrument, such that the 
navigation system determines the position (meaning the location and 
orientation of the medical instrument) relative to the body of the patient. 
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Based on the position determination, thus based on a determination of the 
spatial coordinates of the medical instrument, an image of the medical 
instrument can be bl e nd e d folded into the image associated with the volume 
data set (that is, for example, shown on a monitor). 

An automatic recognition of the approximate position of the medical 
instrument also is us e ful suitable with in nuclear medicine images . A volume 
data set acquired by PET exhibits a relatively good contrast, however it has a 
relatively poor resolution. An alternative is to manually designate the 
geometric center with a mouse click in another representation, for example in 
an x-ray image produced using contrast means in two projections, or in a CT- 
MPR multi-planar reformatting representation. This is nothing more than the 
reassembly of the volume data set in another orientation that, for example, 
originally horizontal slices are re-oriented (orthogonal MPRs: horizontal, 
sagittal, coronal; free MPR: arbitrary slanted secondary slices). This works 
isotropically with new mult i- lin e multi-slice spiral GTs CT or MR devices, for 
example from 512 CT slices with 51 2 2 pixels per slice). The background for 
this is that various views of the image of the subject and its surroundings are 
shown coupled in a defined position and orientation. If the three-dimensional 
view is shifted and/or rotated, the associated MPR representation oriented 
thereto moves with it. The three-dimensional shifting and rotation are 
continued until the image of interest is centrally intersected by the MPR slice. 
Then the middle of this planar representation of the image is clicked; the 
subject center is thus defined in three dimensions. 
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By means of the inventive method, the dimensionality of the volume 
processing (3D image processing) is reduced by at least one stage. The 
curve of a 3D subj e ct object density D of the three-dimensional image along a 
coordinate axis x of the 3D space can be graphically represented as a two- 
dimensional function D1(x). The three-dimensional image in space therefore 
can be correspondingly described with a function D3(x,y,z) whose values 
(characteristic variables) are plotted in the fourth dimension. Thus the space 
spanned by x, y and z can be projected such that only the problem similar to 
the "one-dimensional wand e r i ng random walk " (see Tittmann P., "Irrfahrten in 
Graphen", Memo Hochschule Mittweida, Sept. 1998) or, respectively, an 
optimization along an abscissa x' alone (instead of (x, y, z)) is to be solved. 

In order to achieve this, the planes of the first set of planes are placed 
through the three-dimensional image. In the general case, the planes of the 
first set of planes should [[If]] at a ll possible not least redundantly, but and 
completely, cover all voxels of the three-dimensional image in various 
orientations. The planes of the first set of planes preferably intersect in the 
first straight line that, at least substantially, proceeds through the geometric 
center of the three-dimensional image. With regard to a fixed normal axis, the 
individual planes of the first set of planes can exhibit, for example, equidistant 
angles ai through a m , when the first set of planes comprises m planes. By the 
transition m to a rotat i on rotational geometry, redundancies arise in the 
direction of the rotation point axis , and lessor scanning sparser sampling 
results in the opposite direction. However, scanning primarily needs to be 
observed in the region of the expected surface contour. By the "new 
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acquisition" (re-sampling) of the three-dimensional image in planes rotating 
with respect oppos i t e to one another that intersect the three-dimensional 
image, a first stage of a dimension d e gradation reduction is achieved, since 
the data sets associated with the planes are two-dimensional data sets and, 
for example, are described with Cartesian coordinates. If the volume data set 
is described with Cartesian coordinates (x,y,z), each data set associated with 
the planes of the first set of planes thus can be described with two- 
dimensional Cartesian coordinates (x/'.yj"), whereby j=1 ... m, when the first 
set of planes contains m planes. Subsequently, the data sets associated with 
the individual planes of the first set of planes are analyzed in order to 
determine contours that are imaged in each of the planes of the first set of 
planes and are associated with the surface contour of the three-dimensional 
image. 

In each of these planes (Xj",yj"), a closed contour is expected. The 
contour Kj" of the j-th plane of the first set of planes has the points Kjj" = 
(Xj/'^jj"), with i = 1 ... n. In order to determined each Kjj", a two-dimensional 
search is required in (xj/'.yj/') that is known in general to those skilled in the 
art. 

In order to reduce the calculation effort, however, in an embodiment of 
the invention, when the image contents of each plane of the first set of planes 
are specified with Cartesian coordinates, the following steps are implemented 
for each plane for the determination of the contours of each individual planes. 
A coordinate transformation i s imp le m e nt e d according to polar coordinates is 
implemented with regard to the geometric center of the three-dimensional 



8 



image, and therewith unwinding of the contour, and the contour is determined 
in the transformed plane. 

The contour in the transformed plane is thereby determined according 
to a preferred variant of the invention, in which a contour continuation e nsu e s 
is calculated in the direction of the angular coordinate (O) with an optimization 
according to minimal changes. 

Thus a transformation according to polar coordinates is applied with 
regard to the geometric center and the first straight line. [[A]] Thus a "contour 
field" K'(x',0) with fixed <J> thus results,, that covers 360°. The radial x' is 
designated as r (r=x') in the following, corresponding to the polar coordinates. 
With relatively little effort, only the one-dimensional problem remains, to 
determine the radius r at each angle O as a single quantity. The fact that this 
is actually a rather "one-dimensional" problem can be illustrated by plotting 
this function on a Cartesian system, for example with O in the abscissa and r 
in the ordinate. Some boundary conditions now result for r that can be utilized 
in order to find a qualitatively good contour: 

The contour must have the same r-coordinate for <t>i and O n 
(closed contour). 

The contour is continuous and has a predominantly vertical 
gradient. The r-offset from line to line is relatively small. 
A typical "edge profile" in the brightness profile can be found by 
line comparison. 

Possible compensation of signal no i s e s noise in the line 
comparison. 
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For the transformation according to polar coordinates, the angular 
increment should preferably be selected such that adjacent points (Oj, r) and 
(<tyM , r) in the region of the contour have approximately the samo separation 
distance of a voxel edge length. In place of the transformation accord i ng to 
polar coordinates, a transformation with regard to a "model contour" can also 
onsuo be implemented when a similarity to such a model is to be expected. 
The abscissa is then characterized by the continuous indices tf- ^ through m i^, 
of the contour points of the model; the perpendicular distances to s e parat i ons 
ef the model contour are plotted on the ordinate. Here the a resampling also 
is used, as in the transformation according to polar coordinates is a 
r e scanning . 

In the matrix resulting from the r e scanning resampling , various features 
for an object edge can now be determined or accentuated by suitable signal 
processing, as well as being brought into context. For example, a search 
method can be applied that analyzes and combines all reasonable path 
branchings. An example of a realization is the dynamic optimization 
according to R. Bellman, "Dynamic programming and stochastic control 
processes", Information and Control, 1(3), pages 228-239, September 1958. 
The result is a continuous contour, for example in the polar coordinate 
representation of Oi through O n (or, respectively, for h through i n ). This 
contour can be transformed back Into to the original Cartesian coordinates. 
The dynamic optimization can combine a number of feature values (filter 
responses) into a suitable cost function, and, according to after the forward 
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accumulation of costs, includes a backtracking as a substantial element. A 
plane of the first set of planes is thereby contoured. 

The remaining planes of the first set of planes are correspondingly 
analyzed. Already-determined contours of previously analyzed planes can be 
used in order to also ensure the consistency of the three-dimensional surface 
contour of the three-dimensional image in the rotation direction described by a 
(usable prior knowledge). Likewise, all planes of the first set of planes can be 
analyzed independently. With one or more successive processing stages with 
which the individual contours of the planes of the first set of planes are 
compared, the obtained surface grid can be refined. 

At the end of the processing of all planes of the first set of planes, the 
surface contour of the three-dimensional image initially exists in coordinates 
(r, a, O). In this artificial " cuboid", a processing of the 3D connect i on 
[corr el at i on; i nt e rr e lation] connectivity (correlation, interrelation) is possible, for 
example as a filtering of the r-values with a 2D operator. 

The transformation back ifrte to the original Cartesian coordinates 
yields the surface grid of the three-dimensional image. The surface grid is 
comparable to a net-like cage. The surface grid now can be improved and 
refined in the original coordinates. The voxels within the surface grid can 
subsequently be marked as belonging to the three-dimensional object image. 

The inventive method in the context of a spherical geometry has the 
advantage that it possesses no preferred direction and that, considered 
statistically, in pract i co occurr i ng surface grad ie nts forms that occur in practice 
are handled better on average than in a method that is based on a Cartesian 
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grid, since such a surface lies more parallel to a closed 3D contour. Because 
the inventive method is, geometrically, relatively well-adapted to the problem, 
the number of the planes of the first set of planes can be smaller by a multiple 
than when, in the region of the surface, the surface would be sampled in the 
voxel s e parat i on distance ; this means the "net lines" can have a separation of, 
for example, ten or more voxel e dg e s widths from one another while still 
describing the overall contour with sufficient precision. 

In particular for an improved determination of the surface contouring of 
the three-dimensional image, or when unfavorable output (source) data are 
present, according to an embodiment of the invention the following additional 
steps are implemented. 

A second set of planes is placed through the three-dimensional image, 
with the planes of the second set of planes intersecting in a second straight 
line. The contours that are imaged ar e det e rmined i n e ach of th e p l an e s of jo 
the second set of planes arre determined in each of these planes and are 
associated with the surface contour of the three-dimensional image. Together 
with the contours determined i n oaoh piano of for the first set of planes, the 
contours determined in each plane of the second set of planes are m e rg e d 
combined into the surface grid associated with the surface contour of the 
three-dimensional image. 

By re-extraction oriented on net lines (scan lines proceeding at right 
angles to one another, with orientation from the inside out) determined by 
means of the planes of the first set of planes and associated with the surface 
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grid, the contouring also works with less convex, irregular three-dimensional 
images. 

A benefit ef as to reliability and precision can be achieved when, 
according to an embodiment of the invention, the first straight line is aligned at 
a right angle to the second straight line, and in particular the first straight line 
and the second straight line intersect in the geometric center of the three- 
dimensional image. 

The precise surface contouring can be used, for example, to determine 
the three-dimensional image precisely for a functional evaluation. For 
therapeutic or surgical procedures, the precise shape in relation to the vicinity 
(for example vessels), and also the shape change in the course of time, is 
very important. Further applications are possible in plastic surgery or in 
prosthetics, for the recognition of the subject object to be replaced or the 
recognition of cavities for an implant. 

A simplified calculation of the necessary calculation operations 
illustrates some of the basic advantages of the inventive method. It is 
assumed that the volume data set exists in the form of 50 computed 
tomography slices of the matrix 512x512. The three-dimensional image is 
formed by a sub-cuboid of 40x200x200 voxels. 

Given a surface contouring of the three-dimensional image according 
to the prior art, starting from a start point its continuation on the surface is 
sought. The surface of the entire sub-cuboid is 112000 pixels. For the 
gradient formation, a factor of 1 1 is applied with regard to the surface ( coup le r 
fieW kernel size 3 3 =>"26"), and applied with a factor of "1" for the 27 covered 
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voxels. This effectively corresponds to five surface search-stripes inside and 
outside ef the surface. Starting from each point, in principle 17 directions in a 
preferred direction and laterally are examined (start) that th e n r e duc e ins i d e a 
surfac e initially the number of which then decreases . In a more sensible 
manner, the search is empirically delimited as specified, for example, in 
Buckley K., Vaddiraju A., Perry R., "A New Pruning/Merging Algorithm for 
MHT Multitarget Tracking", Radar-200, May 2000. Furthermore, it is assumed 
that it is sufficient to on average examine six alternatives with an average 
depth of three in order to differentiate a reasonable continuation from a 
m e ander wrong way or, respectively, to land on an already-classified point 
(multiplier 6 3 = 216). 

Given a surface contouring according to the already specified method, 
and not under consideration of "backtrackings", approximately 11x112000x26 
= 33264000 gradient format i ons calculations are required, and approximately 
216x112000 = 24192000 ca l cu l at i on operations for the search for maximum 
gradient sum. Thus approximately 57000000 calculation operations result. 

For the approximate calculation of necessary calculation operations for 
the surface contouring of the same three-dimensional image according to the 
inventive method, it is assumed that the first set of planes comprises 20 
planes. After the transformation into polar coordinates, the image information 
of the transformed planes are respectively sampled in the a 2° raster, 
whereby the corresponding CT images are reconstructed in the typical raster 
512x512: 180x256 = 46000. The gradient formation and contour search are 
first i mpl e ment e d carried out in the planes. Moreover, a noise-suppressing 
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fold i ng (fo l d - ov e r) op e rat i on convolution is implemented with a coup le r fi el d 
kernel size of seven. A typical density profile in the region of the subject 
object edge is generated with seven supporting locations and shifted over the 
entire radius. In each shifting position, the covariance is calculated as an 
i t e rat i on va l u e. Giv e n in the search of the contours, r e sp e ctiv el y a "jump" by 
one pixel from line to line is provid e d allowed ; there are thus three 
alternatives. For post-processing of the connection between the planes, a 
3x3 coupl e r fi e ld kernel is used. It thus results in approximately 20x46000 = 
920000 operations for the gradient formation, and approximately 3x46000 = 
138000 calculation operations for the search for a maximum gradient sum. 
Approximately 20x180x9 = 32400 calculation operations result for the 
backtracking, which here corresponds to a testing of for the plan e plane-to- 
plane connection. Thus approximately 1090400 calculation operations result. 

Based on the inventive method with the c i t e d specified assumptions, a 
reduction of the effort in calculation operations results by approximately the 
factor of 50. 

DESCRIPTION OF THE DRAWINGS 

Figure 1 shows a computed tomography apparatus operable in 
accordance with the invention. 

Figure 2 shows a set of planes whose planes intersect with a three- 
dimensional image acquired with the computed tomography apparatus shown 
in Figure 1. 

Figure 3 shows, in Cartesian coordinates, image information of one of 
the planes shown in Figure 2. 
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Figure 4 shows, in polar coordinates, image information of the plane 
shown in Figure 3. 

Figure 5 is an illustration of a determination of a contour shown in 
Figure 4. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Figure 1 schematically shows a computed tomography apparatus with 
an x-ray source, from which is emitted a pyramid-shaped x-ray radiation beam 

2 (the edge rays of which are shown as dash-dot lines in Figure 1) that 
proceeds through an examination subject (for example a patient 3) and strikes 
a radiation detector 4. In the exemplary embodiment, the x-ray source 1 and 
the radiation detector 4 are arranged opposite one another on an annu l ar a 
gantry 5. With regard to a system axis 6 that runs through the middle point of 
the annular gantry 5, the gantry 5 can be rotated on a mounting device (not 
shown in Figure 1) (arrow a). 

In the exemplary embodiment, the patient 3 lies on a table 7, 
transparent to x-ray radiation, which is movably positioned along the system 
axis 6 by means of a carrier device (likewise not shown in Figure 1) (arrow b). 

The x-ray source 1 and the radiation detector 4 thus form a 
measurement system that, with regard to the system axis 6, is rotatable and 
movable along the system axis 6 relative to the patient 3, such that the patient 

3 can be irradiated at various projection angles and various positions with 
regard to the system axis 6. From the ensuing output signals of the radiation 
detector 4, a data acquisition system 9 thus forms measurement values that 
are supplied to a computer 11 that calculates, by means of method known to 
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those skilled in the art, an image of the patient 3 that can be reproduced on a 
monitor 12 connected with the computer 11. In the exemplary embodiment, 
the data acquisition system 9 is connected with an electrical line 8 that has , 
for example, (in a manner not shown) a slip or wiper ring system, or a wireless 
transmission link, with the radiation detector 4, and with an electrical line 14 
with the computer 1 1 . 

The computed tomography apparatus device shown in Figure 1 can be 
used both for sequential scanning and for spiral scanning. 

In sequential scanning, a slice-by-slice scanning of the patient 3 
ensues. The x-ray radiation source 1 and the radiation detector 4 are thereby 
rotated around the patient 3 with regard to the system axis 6, and the 
measurement system comprising the x-ray source 1 and the radiation detector 
4 acquires a plurality of projections in order to scan a two-dimensional slice of 
the patient 3. A sectional image representing the scanned slice is 
reconstructed from the thereby-acquired measurement values. The patient 3 
is moved along the system axis 6 between the scanning of successive slices. 
This event repeats until all slices of interest are acquired. 

During spiral scanning, the measurement system formed by the x-ray 
source 1 and the radiation detector 4 rotates with regard to the system axis 6, 
and the table 7 continually moves in the direction of the arrow b, meaning the 
measurement system formed by the x-ray system 1 and the radiation detector 
4 continually moves relative to the patient 3 on a spiral track c until the region 
of interest of the patient 3 is completely acquired. A volume data set is 
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thereby generated that, in the case of the present exemplary embodiment, is 
coded according to the DICOM standard typical in medical technology. 

In the exemplary embodiment, the surface of a three-dimensional 
image of a subj e ct object (acquired with the volume data set) to be contoured 
is part of the region of interest of the patient 3. In this case, the subject object 
is a tissue variation of the patient 3. The three-dimensional image 21 is 
schematically shown in Figure 2 and can be viewed with the monitor 12. 

In order to contour the three-dimensional image 21 , a first set of planes 
that (in the exemplary embodiment) contains 20 planes Ei through E 2 o is 
placed through the three-dimensional image 21. For clarity, in Figure 2 only 
the planes E1 through E5 are shown. The planes E1 through E20 intersect in 
a first straight line G that, in the exemplary embodiment, proceeds through the 
geometric center S of the three-dimensional image 21 . The geometric center 
is either manually, approximately marked in the three-dimensional image 21 
by an operator (not shown in Figure 1) by means of a mark (not shown in 
Figure 2) that can be m i x e d faded into the field of view of the monitor 12 and 
that can be moved with a computer mouse 13 connected with the computer 
1 1 , or the computer 1 1 is (in the exemplary embodiment) configured such that 
it automatically determines the approximate geometric center S of the three- 
dimensional image 21, and subsequently automatically determines the planes 
Ei through E 2 o- Moreover, the operator in the exemplary embodiment can 
subject the volume data set to, for example, a median filtering with the 
computer 11. The volume data set is thereby significantly smoothed while 
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preserving edges, in order to improve the determination of the geometric 
center. 

In the exemplary embodiment, the planes E1 through E20 rotat e 
oppos i t e are rotated relative to one another, respectively by an angle Aa = 9°, 
and the first straight line proceeds parallel to the x-axis of the DICOM 
coordinate system already defined for the volume data set. If, for example, 
the three-dimensional image 21 has a quarter of the overall size of the 
isotropic volume data set in the exemplary embodiment, the maximum point 
s e paration distance at the surface of the three-dimensional image 21 thus 
corresponds to approximately 20 pixels. 

In each of the planes E1 through E20, slices of the three-dimensional 
image 21 are imaged. Such a section 31, which in the exemplary 
embodiment is fix e d through defined by the plane E1 , is schematically shown 
in Figure 3. Moreover, the image contents of the individual planes E<| through 
E 2 o are described in Cartesian coordinates. The image contents of the plane 
Ei are described with the Cartesian coordinates x" and y" in the exemplary 
embodiment. 

In order to obtain the surface contour of the three-dimensional image 
21, the slices shown in the planes Ei through E 2 o, thus also the slice 31 of the 
plane Ei, are subsequently analyzed, in that the contours of the individual 
slices are determined. The contour of the slice 31 is provided with the 
reference character 32 in Figure 3. 

For an efficient determination of the contours, in the exemplary 
embodiment the image information of each plane Ei through E 2 o is radially 
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newly arranged, by be i ng transformed with regard to the geometric center S 
according to polar coordinates (r, O). The result of such a coordinate 
transformation is shown as an example in Figure 4, for the transformation of 
the image contents in particular of the slice 31 of the plane E1 shown in 
Figure 3. The transformed contour of the slice 31 is provided with the 
reference character 32' in Figure 4. 

The transformed contour 32' is subsequently determined. A limitation 
path (definition path) is thus sought that corresponds to the contour 32' 
transformed into polar coordinates. For this, a direction-focuses filtering first 
ensues that can even be one-dimensional in the r-direction, and that can 
operate so as to both be sensitive to edges and to suppress noise. A one- 
dimensional search for the limitation path, thus for the transformed contour 
32', now is undertaken using filtering according to the principle of dynamic 
optimization, as, for example, specified in R. Bellman, "Dynamic programming 
and stochastic control processes", Information and Control, 1(3), pages 228- 
239, September 1958. Thus the r-coordinates corresponding to the angular 
coordinates <t> between 0° and 360° are sought, with the boundary condition 
that the r-coordinate for <t> = 0° should be equal to the r-coordinate for O = 
360°. In the exemplary embodiment, these r-coordinates are determined by, 
after a relatively significant smoothing, superimposing the first and last five 
lines, and subsequently the edge positions are determined in the smoothed 
line data. This search is shown as an example in Figure 5, wherein the 
transformed contour is again provided with the reference character 32'. The 
lines 40 in Figure 5 are shown thinner than the transformed contour 32' and 
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represent path alternatives found during the determination that, however, 
were again discarded. 

The contours (transformed into polar coordinates) of the remaining 
planes E 2 through E 20 are similarly determined and subsequently considered 
in context. This involves examining whether the continuity of the brightness or 
of the gradients relevant in this direction is plausible. 

The found contours (transformed into polar coordinates) of each plane 
E1 through E 2 o are subsequently transformed back into Cartesian coordinates, 
and after this are merged into a surface grid. The surface grid thereby 
corresponds to the surface contour of the three-dimensional image 21 . 

In the exemplary embodiment, the surface grid is now improved and 
refined in the original three-dimensional Cartesian coordinates. Using the 
surface grid, the surface contour of the three-dimensional image 21 to be 
determined is interpolated up to the point separation. The voxels within the 
surface grid are subsequently marked as pertaining to the three-dimensional 
image 21. 

The volume data set is shown as a three-dimensional data block B'(r, 
ct>, a) by the placement of the planes through E 2 o through the three- 
dimensional image 21 and via the subsequent coordinate transformation 
according to polar coordinates. The determination of the transformed 
contours (such as, for example, the transformed contour 31') is less 
calculation-intensive for approximately convex subjects. A contiguous 
"mountain surface" is sought via O and a. The "mountain surface" locally 
exhibits relatively minor oscillations in the r-direction. 
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The number of the voxels analyzed for the determination of the 
contours is also reduced with regard to a three-dimensional surface 
contouring according to the prior art. If the volume data set contains, for 
example, 50 CT slices with an approximate dimension of 128 pixels of the 
three-dimensional image whose surface contour should be determined (thus a 
processing field approximately 50x256x256 large), approximately a million 
voxels (20 planes, polar coordinate raster 2°, r max = 256) are analyzed. 
According to the prior art, in contrast to this approximately 3 million voxels are 
analyzed. 

In order to obtain the improved surface grid in the exemplary 
embodiment, in addition to the planes Ei through E 2 o of the first set of planes, 
alternatively at least one further, second set of planes (not shown in the 
Figures for reasons of clarity) can be placed. The planes of the second set of 
planes intersect in a second straight line that likewise proceeds through the 
geometric center S and is aligned at a right angle to the first straight line G. 
The slices shown with these planes are again transformed into polar 
coordinates with regard to the geometric center S. The transformed contours 
that are associated with the surface contour of the three-dimensional image 
21 are subsequently determined again, analogously as above. The 
determined transformed contours are subsequently again transformed back, 
and thereafter are merged into the surface grid. 

It is also an advantage that, via the coordinate transformation into polar 
coordinates, the determination of the transformed contours in the rotation 
direction a (azimuthal) and in directions within the planes (O, orbital) can 
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substantially be three-dimensionally implemented as a two-dimensional 
filtering, namely as a smoothing of r(0,a), instead of in the original 
coordinates. 

The volume data set with which the three-dimensional image 21 is 
produced is, in the exemplary embodiment, produced with a computed 
tomography apparatus. The volume data set alternatively can be produced 
with other medical technology devices, such as in particular with an x-ray 
device, a magnetic resonance device, by PET, an ultrasound device, or by 
means of nuclear medicine. Imaging devices or technical devices outside of 
medical technology also can be used. 

The three-dimensional image 21 is an image of a part of the patient 3, 
however, the surface of a non-human three-dimensional image also can be 
contoured by means of the inventive method. 

The first straight line G need not necessarily proceed through the 
geometric center S of the three-dimensional image 21. The planes Ei through 
E 2 o also need not necessarily intersect in the straight line G. 

Although modifications and changes may be suggested by those 
skilled in the art, it is the intention of the inventor to embody within the patent 
warranted hereon all changes and modifications as reasonably and properly 
come within the scope of his contribution to the art. 
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